Are there plasminos in superconductors? 
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Hot and/or dense, normabconducting systems of relativistic fermions exhibit a particular collective 
excitation, the so-called plasmino. We compute the one-loop selbenergy, the dispersion relation 
and the spectral density for fermions interacting via attractive boson exchange. It is shown that 
plasminos also exist in superconductors. 



PACS numbers: 12.38.Mh, 74.62.Yb 



o 
o 

(N 



o 
o 



X 



I. INTRODUCTION 



In relativistic fermionic systems at high temperature and/or high density, there are two types of fermionic excita- 
tions. Besides the ordinary excitation branches of particles and antiparticles, there are additional collective excitations, 
the so-called plasmino and anti-plasmino (T], Q ■ These branches have opposite chirality compared to the ordinary ex- 
citation branches [5| . They coincide with the normal fermionic branches for vanishing momenta. For large momenta, 
, they approach the lightcone and their spectral strengths vanish exponentially. The plasmino (anti-plasmino) branch 
^ ■ has a minimum (maximum) for smalL non-zero momenta. These excitations have been extensively investigated for 
Q^ ; normal-conducting matter @, i, i, i, 0, i, i, 0, [HI, E [S H [H El, [13, 111 , for a review see, for instance, Ref. [H . 
T— I , With the exception of heavy-ion collisions, in the laboratory it is hard to achieve sufficiently large temperatures 
and densities such that a relativistic description for fermions becomes necessary, even if they are as light as electrons. 
On the other hand, there is a plethora of astrophysical situations where fermions have to be treated relativistically. 
For instance, the core of compact stellar objects could be sufficiently dense to consist of deconfined quark matter. 
The quark Fermi energy is then of the order of ~ 500 MeV. Thus, at least the light up and down quark flavors 
have to be considered as relativistic particles, since m ~ 5 MeV <C fi. However, quark matter in compact stellar 
^ I objects, if sufficiently cold, is not a normal-conducting system, but a color superconductor (20l. [2l|. In this paper, we 
I r. therefore investigate whether the plasminos known from normal-conducting systems survive in a superconductor. This 
' is a question of general interest, independent of the specific nature of the superconductor (i.e., ordinary, color, etc.). 
^ , To our knowledge, this problem has not been considered previously, since electrons can to good approximation be 
C • considered as non-relativistic in the condensed-matter context, and plasminos are absent if the temperature is smaller 
than the mass (at least for zero chemical potential Q). It is also possible to formulate our expectation regarding 
the existence of plasminos in superconductors: plasminos are low- momentum excitations which, for large chemical 
potential fj, ^ p, are buried deep down in the Fermi sea. On the other hand, superconductivity is a Fermi-surface 



■ phenomenon. We thus expect that superconductivity should not exert a destructive influence on the presence of the 
■ " ' plasmino excitations. Nevertheless, it requires an explicit calculation to prove this, which is the purpose of the present 
paper. We shall see that our expectations regarding the existence of plasminos in superconductors are confirmed. 

The outline of this paper is the following. In Sec.|TTl we consider a normal-conducting system consisting of massless 
fermions interacting via scalar and vector boson exchange. For the sake of simplicity, we restrict our consideration to 
zero temperature. This case has been studied before by Blaizot and OUitrault in Ref. [§]. We largely confirm their 
results and extend them by computing the spectral density. We then study a superconducting system in Sec. IIIII 
Section [IVI concludes this paper with a summary of our results. 

Our units are h — c — ks = ^ and the metric tensor is g^'^ = diag(+, — , — , — ). Four-vectors are denoted by capital 
letters, K = (fco, fc). Three- vectors have modulus fc = |fc| and direction k = k/k. Our computations are done in the 
imaginary-time formalism where space-time integrals are denoted as = /q^^^ dr Jy (fix, while energy-momentum 
sums are j^. = Tj^n I d^k/{2Tr)^ with n — 0, ±1 ± 2, ... labeling the Matsubara frequencies for bosons, cu'^ = 2n7rT, 
and fermions, = (2n -I- l)7rT, respectively. 
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II. NORMAL-CONDUCTING FERMIONS 

In this section the dispersion relation and the spectral density of massless normal-conducting fermions is reviewed 
in the limit T @ . We first investigate the case where the interaction between the fermions is mediated by scalar 
bosons and compute the one-loop fermion self-energy. We then use the Dyson-Schwinger equation for the fermion 
propagator to perform a resummation of the one-loop self-energy to all orders. From the thus obtained quasiparticle 
propagator we determine the dispersion relation and the spectral density. Finally, we conclude this section with a 
discussion of the case where the interaction is mediated by vector bosons. Note that this treatment is not fully 
self-consistent in the sense that we do not use the result for the resummed quasiparticle propagator to recompute the 
one-loop self-energy for an iteration of the above procedure. Moreover, we do not solve a Dyson-Schwinger equation 
for the scalar boson propagator, which would lead to the generation of a boson mass ~ gfi. 

A. Self-energy 

The inverse propagator for massless non-interacting fermions can be written in the form 

Go i(P) = f + ^,Jo ^ 70 E ^o.liP) , (1) 

e=± 

where we introduced the energy projectors 

A± = i(l±707-j5), (2) 

and the free inverse propagators for positive/negative-energy solutions 

GoX{P)^Po + ^^TP. (3) 

Equation Q can be easily inverted to give 

Go(P) = EGo,e(P)Al7o, (4) 

e=± 

with Go,±(P) = l/(po + ^iTp)• 
For fermions interacting with scalar bosons, the interaction part of the Lagrangian is of Yukawa-type, Ci = gipij^cj). 
The (perturbative) one-loop self-energy has the form 22] 

nP)^-9^TY, -^Vo{K-P)Go{K). (5) 

We assume the boson to be massless as well, with propagator 'Dq{Q) = —1/Q^, thus there is no tadpole contribution 
to the fermion self-energy. Performing the sum over the Matsubara frequencies one obtains 



l-N^p{k)+NBiEb) , N^p{k) + Nb{Ei,) 



po + - e{k + Eb) po + - e{k - Eb) 



(6) 



Here, Eb = \k ~ p\ is the energy of the exchanged boson, Np{E) = [e^^^'^^/^ -t- are the thermal distribution 
functions for fermions and antifermions, and Nb{E) = [e^l'^ — is the corresponding one for bosons. 
We now project onto the self-energies for positive/negative-energy solutions, 

E±(F)^iTr[A±7oS(/')] , (7) 

and perform an analytic continuation, po + /i ^ a; + 177, where io is the (real-valued) fermion energy relative to the 
vacuum 

Using Eq. ([7]), one readily calculates the imaginary part of the self-energy for positive/negative-energy solutions, 



ImE±(a;,p) — Ima(w,p) ±plm6(a;,p) 



(8) 



FIG. 1: The imaginary parts of the self-energy for massless fermions (left panel) and antifermions (right panel) at T = and 
for an exemplary value of the coupling constant /{An) — 1. 

I Re(atb'p)/n | | Re(a-b•p)/^ | 




FIG. 2: The real parts of the self-energy for massless fermions (left panel) and antifermions (right panel) at T = and for an 
exemplary value of the coupling constant /{Air) = 1. 

Due to rotational invariance, the functions lma{u!,p) and Im6(w,p), which were first calculated in Ref. only 
depend on the modulus of the momentum. For the sake of completeness we list the results in Appendix \^ The real 
part of the self-energy for positive/negative-energy solutions is 

Rel]±(ci;,p) — Rea(w,p) ±pRe6(a;,p) , (9) 

with the functions Rea(cij,p) and Re6(a;,p) given in Appendix \K\ 

We show the functions ([5]) and ^ in Figs. [1] and [21 respectively. The shape of the imaginary and real parts reflect 
the different regions shown in Fig. 1131 The peak of the imaginary parts in the region of low energies and momenta is 
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due to the 1/p singularity of the function Ima, cf. Eq. (jA5p . The corresponding peak seen in the real parts, however, 
is due to the singularity of the function Re 6, cf. Eq. (|A11[) . This is the reason why the singularities appear with 
the same signs in the imaginary parts, but with opposite signs in the real parts of the self-energies for particles and 
antiparticles, cf. Eq. ([7]). 

B. Dispersion relation 

The full inverse fermion propagator is given by 

G-\P) = Go\P) + ^ 70 E G-HP)AI, , (10) 

e=± 

where the full inverse propagator for positive/negative-energy solutions is 

G^\P)^G,XiP) + ^±iP) ■ (11) 
The dispersion relations po + ^ = lu = i-u'^{p) are given by the roots of the real parts of the inverse propagators, 

ReGz\Lol,p) = 0. (12) 

We show the solutions of these equations in Fig. [3l The ordinary particle and antiparticle excitations correspond to 
the uppermost and lowermost curve, respectively. Besides these, for low momenta we find two additional roots of 
both Eqs. (fT2|) . The two solutions in the time- like region belong to the plasmino and antiplasmino branches. For 
massless particles right(left)-handedness implies positive (negative) helicity, while for massless antiparticles right(left)- 
handedness implies negative (positive) helicity. Since plasminos have opposite chirality from particles, the plasmino 
solution is actually a root of GZ^ — 0, while the antiplasmino is a root of G^^ = 0. This also holds for the two 
additional solutions in the space-like region. However, since the imaginary parts for particles and antiparticles are 
large in this region, these excitations are strongly damped and do not have appreciable spectral weight, cf. Fig. 3] 
The plasmino solutions approach the lightcone and can no longer be found numerically for large momenta. The size 
of the momentum region where plasminos are found depends on the value of the coupling constant. For decreasing 
coupling constant, this region shrinks. 



C. Spectral density 



The spectral density is determined from the relation [25 



1 



Im G±{uj,p) 



(13) 



If the imaginary part of the self-energy vanishes (see region IV of Fig. I13p . the spectral density is proportional to a 
^-function with support on the quasiparticle mass shell, u = LUj.{p), 



P±(w,p) =E ^±(wi)(^(cJ-cj±) , Z±{uj^) 



dRc G±\uj,p) 



(14) 



which corresponds to an infinite lifetime of the associated quasiparticle. This is the case for the plasmino branch and 
that part of the particle excitation branch, that is shown in Fig. [31 (For larger momenta, the particle excitation branch 
enters region la of Fig. [13] where the imaginary part is non-zero and, consequently, the particle excitation becomes 
unstable.) On the other hand, a non-zero imaginary part gives rise to a non-zero width of the spectral density. 



P±(w,p) 



lTaT,±{uj,p) 



^ [Go-^(a;,p) + ReI]±(^,p)] + [Iml]±(a;,p)]^ 



(15) 



This leads to a finite quasiparticle lifetime which is inversely proportional to the width of the spectral density around 
the quasiparticle mass shell. This is the case for all other fermionic excitation branches which lie outside region IV of 
Fig. [131 
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FIG. 3: Dispersion relations of fermions for g /(iir) = 1 at zero temperature. The thin lines represent the lightcone. 
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FIG. 4: Contour plot of the spectral density of fermionic excitations in the limit of vanishing temperature and for g /(47r) — 1. 



In Fig. |4]we show the spectral density in a contour plot. Comparing this figure to Fig.[3l one can easily distinguish 
the particle, plasmino, antiplasmino, and antiparticle branches (from top to bottom). As mentioned above, the two 
solutions in the space- like region do not have sufficiently large spectral weight to appear in the contour plot. In other 
words, the width of the spectral density at the respective mass shell is large and, consequently, these excitations 
are strongly damped. Note also that the spectral weight of the two plasmino branches rapidly decreases for larger 
momenta. 
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D. Fermions interacting via vector bosons 



The self-energy for fermions interacting via vector boson exchange is 

^V,^{K-P)rGo{K)j'^ . (16) 

For massless vector bosons, this is a gauge-dependent expression. One easily convinces oneself that, in Fcynman gauge 
for the gauge boson propagator, Sv(^') — 2 S(P). It was shown in Ref. Q that the same result holds also in Coulomb 
gauge. Finally, in the hard-thermal-loop (HTL) and hard-dense-loop (HDL) limit i.e., for T,fi^ Po,P, the result 
is gauge-independent, and again we have Sv(^') = 2I](P), provided one takes the HTL/HDL limit on both sides of 
this equation. Therefore, all these cases are covered by the above discussion for scalar boson exchange via rescaling 
the coupling constant, g — > g/\/2- 



III. SUPERCONDUCTING FERMIONS 



We now turn to the question whether there are plasminos in superconductors. Our rationale is the following. We 
take the superconducting ground state as a basis for our calculation. Then, the fermion in the loop in Eq. ([5]) is no 
longer represented by a free propagator, Go{P), but by a quasiparticle propagator appropriate for superconducting 
systems, Gq. We compute the one-loop self-energy for quasiparticle excitations in the superconductor. Afterwards, 
we resum the one-loop self-energy via a Dyson- Schwinger equation to obtain the propagator for quasiparticles and 
charge-conjugate quasiparticles. The poles of the propagator determine the dispersion relation. Finally, we calculate 
the spectral density. 

Note that this treatment is not fully self-consistent. As in the normal-conducting case, we do not use the resummed 
quasiparticle propagator to obtain a self-consistent result for the one-loop self-energy. Moreover, we do not solve a 
gap equation for the superconducting gap parameter. Rather, we shall use an exemplary constant value for the gap 
function of particles, 4>+{P) — 0.25/^ and set the gap function for antiparticles to zero, 4>-{P) ~ 0. The assumption 
of a constant gap function is justified only for point-like four-fermion interactions as in e.g. the Nambu-Jona-Lasinio 
model. For a boson exchange interaction of non-zero range, the gap function does depend on energy and momentum. 
Setting the antiparticle gap to zero is justified only in weak coupling, because then the contribution of the antiparticle 
modes to the gap equation is negligible [2^, '25']. Since antiparticle modes live far from the Fermi surface, a non-zero 
value (of reasonable magnitude) for the antiparticle gap would only have a negligible influence on our results for the 
excitation spectrum and spectral density of antiparticles. 

Finally, we also do not solve a Dyson-Schwinger equation for the scalar boson propagator which would give rise to a 
mass ~ gp,. Nevertheless, we expect that the complications arising from a fully self-consistent treatment, including the 
energy-momentum dependence of the gap function and non- vanishing antiparticle gaps, will only lead to quantitative, 
but not qualitative, changes of our results concerning the existence of plasminos in superconductors. 

Let us finally note that, in a superconductor, fermions carry a charge with respect to a local (gauge) symmetry. 
This symmetry is spontaneously broken by condensation of Cooper pairs. An implicit assumption of our approach is 
that interactions due to exchange of vector bosons of the gauge interaction are negligible as compared to scalar boson 
exchange, i.e., that the one- loop fermion self-energy is dominated by the scalar interaction. If the fermions do not 
carry any charge, we are strictly speaking not considering a superconductor but a superfluid. 



A. Self-energy 



If fermions interact via scalar boson exchange, and if we only allow for pairing in the even-parity channel, the 
quasiparticle propagator reads [1^ [2^ 

(Note that our Qq corresponds to in Ref. p5l|.) 

The analogue of Eq. ([5]) for the superconducting system now reads 

S(P) = / ^'^oiK-P)GoiK). (18) 
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We compute the Matsubara sum using the mixed representation for the boson propagator, 

V^ir.q) ^ ^^{[1 + NB{E,)]e-^^^ + NB{E,)e''-^} , (19) 
as well as the quasifermion propagator, 

e=± 

where te = \/ (/-* — e^)^ + |0ep, n'^p = [e'^^/"^ + At vanishing temperature, the self-energy then becomes 

We now project onto positive/negative-energy states according to Eq. ([7]). After analytic continuation, the results for 
the imaginary and the real part of the self-energy can again be written in the form of Eqs. ([5]), The functions 
Ima(ci;,p), Im6(cj,p) are available in closed form and are given in Appendix [Bl However, the real parts Rea(w,p), 
Re6(a;,p) have to be computed numerically. The results are shown in Figs. [5] and [S] 

Again, the shape of the imaginary and real parts reflects the different energy domains in Fig. [Ml which are slightly 
different from those in Fig. 1131 Apart from this, the overall shape is rather similar to that of the corresponding 
imaginary and real parts for normal-conducting fermions, cf. Figs. [1] and [2] Note, though, that the peaks are somewhat 
broader and flatter. This means that the damping of quasiparticles is larger than in the normal-conducting case. 



(1 



2e, 



+ {n — ek) 
2Z 



A|7o 



(20) 



/ 1 ^ A^M^A 1_ 

V ee J Po- f^e- Eb \ ee J Pa + Se 



B. Dispersion relation 



We now compute the dispersion relation for superconducting fermions. As usual, it is advantageous to distinguish 
particles and charge-conjugate particles by introducing the Nambu-Gor'kov basis. The propagator for quasiparticles 
and charge-conjugate quasiparticles then reads, cf. Eq. (136) of Ref. plj . 

e± = (^[Gj]-i + S±-$T{[Gi[]-i + ET}"'$±^ . (22) 
Here, the inverse free propagator for particles and charge-conjugate particles is 

[Gt]-\P) =loJ2 N ± - ■ (23) 

e=± 

Note that is identical with Go of Sec. |TT1 Also, if we set = 0, Q'^ becomes identical with in Eq. (|17p . 
Furthermore, the self-energy for particles, S"*", is the one we computed above, cf. Eq. ^TE\i . In the following, wc only 
need its decomposition in terms of projections onto positive and negative energies, 

S+(P)-7oE^^(^)^^ (24) 

e=± 

The self-energy for charge-conjugate particles, S", can then be computed from S+ via [2l| 

S-(P) ^ C[S+(-P)]^C-i = -70 E Se(-P) A- , (25) 

e=± 

where C = 17^70 is the charge-conjugation matrix, and where we have used the identity C^^C~^ = — 7J. Finally, 
<f>^ is the order parameter for condensation. If condensation occurs in the even-parity channel [25l |. 

'S>+{P)^Y.MP)^pl5- (26) 
The charge-conjugate order parameter <I>~ is defined as 

(P) ^ 70 (P)] So - - E ['^^ (^)] * 75 . (27) 

e=± 
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FIG. 5: The imaginary parts of the self-energy for massless superconducting fermions (left panel) and antifermions (right 
panel), at T = 0, for <j}-t- = 0.25 /i and </!>_ — 0. 




FIG. 6: The real parts of the self-energy for massless superconducting fermions (left panel) and antifermions (right panel), at 
T = 0, for 0+ = 0.25 n and (f)- = 0. 



The dispersion relation is given by the poles of the propagator ([221) • Inserting Eqs. 
obtain 



((7f|) into Eq. 1^, we 



e=± 



(28) 



In order to see which poles belong to solutions for positive or negative energy, respectively, we perform the following 



0.5 1 1.5 2 0.5 1 1.5 2 

P/H p/|i 



particle hole ■ plasmino plasmino-hole 

antiplasmino antlplasmino-hole antiparticle antihole 

add. coll. exc. add. coll. exc. 



FIG. 7: Dispersion relations for superconduct- 
ing fermions with positive energy for T — 0, 
g^/(47r) = 1. The thin lines are the excitation 
branches for quasiparticles and quasiparticle-holes 
when setting = 0. 



FIG. 8: Dispersion relation for superconduct- 
ing fermions with negative energy for T — 0, 
(j^/(47r) = 1. The thin lines are the excitation 
branches of normal antiparticles and antiparticle- 
holes. 



projection of the quasiparticle propagator 



t;±(F)^-Tr[t;±(P)7oA;-] 



This gives 
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QZ{P) 











^P- 


s+ 


i-P) 








bo - 




-P- 
1 


^S+(P)][po 


- 


f p 




-P)] 






Po -\ 




1 


S_(P) ' 














Po - 


- /i - 


-P- 


Po+ 


-p- 


f s 


^P) 








[Po - 


f M 


-P- 


^S+(P)][po 


- fi- 






~P)] 




_(P)P 



(29) 

(30a) 
(30b) 
(30c) 
(30d) 



Since we have set (j)- — 0, there is a cancellation of terms between numerator and denominator in the propagators 
and Q^. This leads to a reduction in the number of poles. 

After analytic continuation po + ^ u + irj, we obtain from the poles of Q'\_ and QZ (these propagators have the 
same poles) the dispersion relation for particle, hole, antiplasmino, and antiplasmino-hole excitations (remember that 
the latter belong to the positive-energy part of the spectrum, and remember also that charge-conjugate aTifi-particles 
correspond to positive-energy solutions). These dispersion relations are shown in Fig. [7] for a value of the coupling 
constant / [Air) = 1. For comparison, the thin lines in this figure show the excitation branches of quasiparticles when 
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setting E+ = 0, i.e., the poles of the propagator (fT7|) for the upper choice of signs. As seen in the normal-conducting 
case, Fig. [31 there were additional collective excitations in the space-like domain. We also find the corresponding 
excitations (plus their hole counterparts) in the superconducting case. These excitations do not have an appreciable 
spectral density (see below), they are strongly damped. The gap in the excitation spectrum is clearly visible. However, 
compared to the gap 2 = 0.5 fi in the case = (thin lines), it is shifted with respect to the Fermi surface p — 
to smaller values of momentum. This effective shift of the Fermi surface is well-known from Fermi-liquid theory and 
decreases with the coupling constant g. 

From the poles of we obtain the dispersion relation of antiparticle and plasmino excitations (as well as that of 
the additional excitation) , while the poles of correspond to plasmino- hole and antiparticle- hole excitations (as well 
as to the additional excitation in the hole sector) . These dispersion relations are shown in Fig. [8l For comparison, 
we show the excitation branches of normal antiparticles and antiparticle-holes. As in the normal-conducting case, 
plasmino excitations as well as the additional strongly damped ones exist only in a region of momenta smaller than 
the Fermi momentum. The size of this region shrinks with decreasing values of the coupling constant. 

C. Spectral density 

The spectral densities corresponding to the propagators pO|) are computed via 

pt{u,p)^--lmgf{io,p), (31) 

TT 

e = ±. For the explicit calculation, see Appendix [Cl 

In Fig. [HI we show a contour plot of pj. For small momenta, the particle and hole branches are undamped, since in 
this region the imaginary part of the self-energy vanishes. Consequently, the spectral density is a delta function. For 
larger momenta, due to a non-zero imaginary part of the spectral density broadens on both the particle and the 
hole branch, however, for the particle branch this broadening is less pronounced than for the hole branch. One also 
observes the antiplasmino branch which is strongly damped by a non-zero imaginary part of the self-energy. Although 
the antiplasmino-hole branch corresponds to a pole of Q^, it does not have appreciable spectral strength in p'^, due 
to a numerical cancellation of terms in the numerator of the expression for p^, cf. Eq. (jCl[) . 

In Fig. [TO] we show the spectral density pZ ■ The spectral strength on the particle and hole branch is comparable 
to that for pj. However, in contrast to that spectral density, now the antiplasmino-hole branch is clearly visible and 
the antiplasmino branch is suppressed due to a numerical cancellation of terms, cf. Eq. (|C4p . 

In Figs . [TT] and [T2l we show the spectral densities pi and p^, respectively. We clearly distinguish the plasmino and 
antiparticle branches, as well as the corresponding hole excitations. The antiparticle and antiparticle-hole branches 
are strongly damped because of the non-vanishing imaginary part of the self-energy in this region of the energy- 
momentum plane. Note that results similar to ours have been obtained in Ref. [27|, applying the Maximum Entropy 
Method (MEM) to extract the spectral functions from Dyson-Schwinger equations. 

IV. SUMMARY 

In this work, we have studied the fermionic excitation spectrum in normal- and superconducting systems. For the 
sake of simplicity, we have taken the fermions to be massless and we have considered only the most simple case of an 
interaction mediated by massless scalar bosons. Moreover, we restricted our consideration to zero temperature. In 
normal-conducting, relativistic fermionic systems, there exist additional collective excitations, the so-called plasminos. 
These have opposite chirality compared to the usual particle excitations. 

The goal of the present work was to investigate whether such excitations also exist in superconducting systems. 
In order to answer this question, we have first computed the one-loop fermion self-energy. Using this self-energy, we 
have determined the resulting dispersion relation from the poles of the full fermion propagator. We do indeed find 
plasmino excitations in superconducting systems. 

In addition, we have also identified collective excitations which, at least for normal-conducting systems, lie in the 
space-like region of the energy-momentum plane. Both the plasmino and these additional collective excitations exist 
only in a region of low energies and momenta. We then calculated the spectral density and found that the additional 
collective excitations are strongly damped and do not carry appreciable spectral weight. 
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FIG. 9: The spectral density for superconduct- FIG. 10: The spectral density p_ for supercon- 

ing fermions, T = 0, g^/(47r) = 1, = 0.25 /.i, ducting fermions, T = 0, g^/i^Tv) = !,(/)+ = 

and = 0. 0.25 /i, and (f>- — 0. 




FIG. 11: The spectral density p_ for supercon- 
ducting fermions, T — 0, (7^/(47r) = 1, (^+ = 
0.25 ^i, and <j>- = 0. 



FIG. 12: The spectral density for supercon- 
ducting fermions, T = 0, (/^/(Itt) = 1, </>+ = 
0.25 /i, and (f)- = 0. 
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APPENDIX A: IMAGINARY AND REAL PART OF THE SELF-ENERGY OF NORMAL-CONDUCTING 

FERMIONS 

The imaginary and real part of the self-energy, see Eqs. ((8]) and ([9]), respectively, can be computed analytically. 
One distinguishes different regions in the energy-momentum plane as shown in Fig. 1131 In these regions, the functions 
Ima(a;,p) and lmb{uj,p) assume the following values, see Ref. [Oj: 




FIG. 13: The different domains in the energy-momentum plane, which occur in the computation of the imaginary part of the 
self-energy. 
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The real part of the self-energy is the sum of a vacuum and a matter contribution. After renormalization one 
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obtains at the renormalization scale /i the result: 



Iiea{uj,p} 



J {2^l - LU - p)i2ii + u; + p) 
{2fi — bj + p) {2fi + uj — p) 



Ap 



In 



In 



U! +p 



cij + p — 2/i 
a; — ]3 



UJ — p — 2^ 



UJ In 



uj^ — p^ 



5^ J 



{2^.-uj-p) 
(2^ — w + p) 



4p3 



uj"" -p" - -{2n + ^+p) 



uj"^ -p^ - -(2/i + LJ-p) 



In 



In 



UJ +p 



cj + p — 2/x 

w — p 



— p — 2/x 



1 



In 



(AlO) 



(All) 



Comparing this to Eq. (B6) of Ref. @, we obtain an additional term —ujp,/{2p'^) in Eq. (jAlip . 
APPENDIX B: IMAGINARY PART OF THE SELF-ENERGY OF SUPERCONDUCTING FERMIONS 



Superconductivity is a consequence of the formation of Cooper pairs at the Fermi surface, i.e., it involves particles, 
but not antiparticles. Therefore, it is permissible to set the antiparticle gap function to zero, 0_(P) = 0. In the 
following, we abbreviate the particle gap function as 4>+{P) = 4>- The fact that particles have a gap but antiparticles 
do not forces us to compute particle and antiparticle contributions to the imaginary part of the self-energy (j2ip 
separately. We therefore split the latter into two parts. Using the projection ([7]), we obtain 



e=± 



Im4^'(^,p) 
Ima'-'^''(a;,p) 

plmb^'^\uj,p) 



lma''^\uj,p) ±plv[\b^'^\uj,p) 
2 I (fk 1 



ng 



(27r)3 8Eb 



1 _ [L^] s{p„ 1 + ) Sipo + 6e + Et) 



eng 



d?k p ■ k 
(2^)3 S-Bfc 



1 - S{po -ee^Eb)+(l + ) S{po + ee + E,) 



(Bl) 
(B2) 
(B3) 

(B4) 



As in the case of normal-conducting fermions, we encounter different regions in the energy- momentum plane when 
computing the functions Ima*^'^^(a;,p), lmb'^'^^{uj,p). For the negative-energy contribution (e = — ) to the self-energy 
(|2T|) . these regions look the same as in Fig. [131 because we have set the antiparticle gap to zero. However, for the 
positive-energy contribution (e = they are of slightly different shape due to the non-zero value for the particle 
gap (f). We show these regions in Fig. [141 

In the following, we list the functions Im a'^^'^^ (a;,p), lmb'^^\uj,p) in the different domains of Fig. [HI 
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lb: ImaW = 
la: Ima^^' 

II: Ima(+) = 

Ini6(+) = 

III: Inia(+^ = 



9 



1 in 



327r 
327r4p 
327r 2p' 



1 



2ti)p \ cij — p — 2/i 

Lu — 2n f Lo + p ~ 2n 
1 m ' 



2p 



a; — p — 2/i 



(2^ - cj - p) (2/^ + + p) + 



4m , r 

uj — p — p)^ 



2 In 



{lj - p){iL> + p — 2fi) 



LU^ - p'' - - {2ii + u; + p) 



UJ ^ p 



'2fip 



02 - 2p 
2 UJ — p 



5' 1 



- 2/i) In 
4/1 , 02 



UJ — p {uj — pY 



2 In 



|^^{(2m---p)[-^-P^-|(2/. + .+p) 

^2 - 2p 



IV: Inia(+) = Im6(+) = 



^^02 
327r^ 

327rp2 



— p 

2/1 

p2 — cj2 



-2/ip 



2 UJ ~ p 



+0^ (lj - 2//) In 



(p-tj)(2/i-Lj-p) 



^In^"-^ 



(p2 — ti;2)2 2p \UJ + p 



02p 



2flUJ 

Up' — p2 (tj2 _ p2'jS 



2uL — UJ / UJ + p 
■ m ' 



2p 



UJ — p 



(B5) 
(B6) 
} , (B7) 

(B8) 
, (B9) 

(BIO) 

(Bll) 

(B12) 
(B13) 



The negative-energy contribution to the self-energy vanishes except in region III of Fig. [T31 where we find 

^2 



III: Ima^-^ 



32tt 



9 
327r' 

2 



(B14) 
(B15) 
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APPENDIX C: SPECTRAL DENSITY OF SUPERCONDUCTING FERMIONS 

The spectral densities for superconducting fermions are proportional to 5-functions in the region of the energy- 
momentum plane, where the imaginary part of the self-energy is zero. In this case, they can be computed from an 
equation analogous to Eq. (I14p . On the other hand, if the imaginary part is non-zero, they can be computed via 

pI{uj,p) = ^ ^2 I ^2 { - 2^ + p - ReS+(^ - lu, ~p)] B + ImS+ (/i - uj, a} , (CI) 

pZ(uJ,p) = 5 o , (C2) 

T^[uj+p + RcJ:^{Lo-fi,p)f + [lm.i:_{u;-p,p)f 

N 1 Iml]_(/i- 

pAuj,p) ^ ^, (C3) 

+ 7r[cj-2^-p-ReI]_(/x-cj,-p)]' + [ImS_(^-c^,-p)]' 

= ^-^j^-^{[tJ-p + ReI]+(w-M,p)] B-ImS+(w-M,p)v4}, (C4) 

where we have abbreviated 

A = (a; — /i)^ — (/i — p)^ + (ti> — 2/i + j3)ReI]+(ti' — + (p — Li;)ReI]-|_(/i ~ a;, — p) 

— ReE+(u; — /i,p)ReE-|_(/i — w, — p) -I- ImE+(ti> — /i, j5)ImS+(/z — w, — p) — (fP' , (C5) 

i? = (tj — 2/i -|-p)lmS+(a; — /i.p) + (p — a;)lmE-|_(/i — cij, — p) 

— Iml]+(w — /x,p )ReI]+(/x — w, — p) — ReE+(ti> — /i,p)ImI]+(/i — — p). (C6) 
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